Statistical properties of swarms of self-propelled particles with repulsions across the 

order-disorder transition 
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We study dynamic self-organisation and order-disorder transitions in a two-dimensional system 
of self-propelled particles. Our model is a variation of the Vicsek model, where particles align the 
motion to their neighbours but repel each other at short distances. We use computer simulations 
to measure the orientational order parameter for particle velocities as a function of intensity of 
internal noise or particle density. We show that in addition to the transition to an ordered state 
on increasing the particle density, as reported previously, there exists a transition into a disordered 
phase at the higher densities, which can be attributed to the destructive action of the repulsions. 
We demonstrate that the transition into the ordered phase is accompanied by the onset of algebraic 
behaviour of the two-point velocity correlation function and by a non-monotonous variation of the 
velocity relaxation time. The critical exponent for the decay of the velocity correlation function in 
the ordered phase depends on particle concentration at low densities but assumes a universal value 
in more dense systems. 



I. INTRODUCTION 



Self-organisation and non-linear dynamics in systems 
of active objects have been receiving much attention from 
the scientific community during the last years, in particu- 
lar in relation to a number of spectacular biological prob- 
lems. One of the most prominent biological examples is 
the swarming behaviour - spontaneous onset of organised 
motion of large amount of independent individuals. Such 
collective motions are commonly observed [l| in schools 
of fish : bacterial colonies Q , slime moulds [1, 0] and 
flocks of birds [1]. Vortex motions were specifically ob- 
served long ago in colonies of ants @ and more recently 
in bacterial colonies [l^l and slime moulds 11 1 . The col- 
lective motion has been shown to arise from simple rules 
that are followed by individuals and does not involve any 
central coordination, external field or constraint. Re- 
cently, enormous attention has been attracted also to 
swarm robotics, a promising new approach to coordi- 
nate the behaviour of large number of primitive robots 
in decentralised manner through their local interactions. 
A novel area of research. Swarm Intelligence [l^l , covers 
the general features of swarming and collective behaviour 
and expands very rapidly. 

Observation of dynamic patterns in ensembles of bi- 
ological and physical species, which can be generically 
referred to as active species, stimulated a significant ef- 
fort in physics of dynamical systems (for which we refer 
the reader to comprehensive reviews [Tsl - flGj ). A min- 
imal model of swarming was suggested by Vicsek and 
coworkers in 1995 |17| . A comparison of subsequent min- 
imal models has been conducted in Ref. UM. In the 



Vicsek model, the swarming results from the action of a 
collision-type interaction that aligns the velocity of each 
actively moving particle in a group to the mean local ve- 
locity. The statistical properties of the Vicsek model in 
various implementations, the pattern formation and the 



type of order-disorder transition have been extensively 
studied by the community [19l - l24|. T he model was fur- 
ther developed by Couzin et al. [25[ and later was en- 
riched with biologically relevant features, including new 
types of behavioural interaction, which allowed one to 
reproduce several types of swarming dynamics [2l|, [2^ . 

An analogy of the velocity alignment in the Vicsek's 
to the alignment of magnetic dipoles interactin g vi a lo- 
cal mean field was explored by Toner and Tu |27l . [28| , 
who predicted an existence of the aligned (ferromagnetic) 
and isotropic (paramagnetic) steady states in a generic 
continuum model of an active system with aligning in- 
teractions. The transition from the disordered behaviour 
to aligned collective motion occurs in all these models 
on decreasing the noise strength, similar to the magnetic 
transition on decreasing temperature. The correspond- 
ing phase diagram, however, is more complex than the 
typical magnetic one as active particles are not fixed on 
the lattice and their concentration can vary across the 
system. The aligning interactions in this case act as a 



cohesive factor [lj|. As a result, the ordering occurs at 



certain minimum particle number density and the active 
particles demonstrate a tendency to form compact or- 
dered flocks. The orientational transition in the Vicsek 
model with aiigular noise was shown to have continuous 
character (23l[2l|. 

As the transition is associated with orientational sym- 
metry breaking in all model systems, one can expect 
that they belong to the same universality class and 
demonstrate the same critical properties across the tran- 
sition [l3|. A large variety of statistical properties have 
been analysed over the recent years and the validity of 
the scaling relations between critical exponents, known 
from equilibrium phase transitions, has been confirmed 
j23l [29I [30j . The universal properties have been found 
in a variety of active systems including those with only 
pairwise interactions [3l[. Beside the critical behaviour, 
other statistical characteristics can be associated with 
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Figure 1: The interaction parameters in the Vicsek-type 
model. The particle turns away from the nearest neighbours 
within the zone of repulsion (the inner circle, zor) to avoid 
collisions and aligns itself with the neighbours within the zone 
of alignment (the outer circle, zoa). 



transition such as cluster size distributions, showing a 
power law decay (isl [20j , and giant number fluctuations 
HH . We should note that some of these features have 
been also observed in active systems without aligning in- 
teractions but with just hard-core repulsions, so they arc 
not unique to the orientational transition [33 . [33j . Some 
other statistical properties like the two-point velocity cor- 
relation function, however, have not been studied in de- 
tail. Therefore, it would be interesting to test whether 
they also demonstrate the universal behaviour, which is 
not sensitive to the details of the particles and their in- 
teractions. 

In this work, we use a Vicsck-typc model with added 
repulsive interaction (collision avoidance) to study its 
complete phase diagram in terms of particle density and 
noise amplitude. We analyse the behaviour of the orien- 
tational order parameter as well as the most important 
correlation functions across the transition and attempt 
to identify the statistical indicators of the orientational 
ordering. In Section 2 we describe the model and statis- 
tical methods, in Section 3 we describe and discuss the 
results, and conclude the discussion in Section 4. 



II. THE MODEL AND METHODS 

A. Model and simulation settings 

In our two-dimensional model N point particles move 
with a constant speed vq inside a square simulation box 
of size L X L with periodic boundary conditions, so that 
the particle number density is given hy p = N/L'^. The 
direction of motion of each particle is affected by repul- 
sive or aligning interactions with other particles located 
at distance equal or smaller than the radius of the zone 
of repulsion (zor ) or zone of alignment (zoa) respectively 
(see Fig. [T]) [25|, [sj. The first behavioural zone, repre- 
sented in Fig. [T]by a circle of smaller radius _Rr, is re- 



sponsible for the maintenance of a minimum distance be- 
tween neighbouring particles, while the second one with a 
greater radius Ra is used for aligning the particle veloci- 
ties. The size of the repulsion zone is always set to unity, 
so it defines the unit of length for the simulation. We 
note that Rr can also be considered as the body cross- 
section diameter for the particle as the penetration into 
it is penalised. 

At each time step, position of each particle (rj) is com- 
pared to the location of the nearest neighbours. If other 
particles are present within the zone of repulsion, the in- 
dividual attempts to keep its personal space and to avoid 
collision by moving away from the neighbours. This rule 
has been proposed by Couzin et al. ([25l. Isl]) 



V,{t + At)^R{Ut))Vl{t + At) 



with 



V[(t + At) 



(1) 



(2) 



where r,- 



is a number of particles inside zor 



and R{^i{t)) is the rotation matrix 



Rim) 



_ ( cos^i(t) -sin^j(t) 
sin^i(t) cos^i(t) 



(3) 



^i(t) is a random variable uniformly distributed around 
[— TyTT, ?77r], where rj represents the noise strength. 

Repulsions have the highest priority in our model, 
therefore no alignment rule can be applied if other parti- 
cles are present in the zor. Only if the zone of repulsion 
contains no neighbours, the particle responds to the next 
zone and to align to other particles' motion 



V,(t + Ai) = i?te(t))VJ'(t + Ai) 



with 



EV,(t) 



(4) 



(5) 



where is a number of particles inside zoa. 

The above-described behavioural rules are used 
throughout the paper, however, in the last section we 
also show (where especially noted) some results for mod- 
ified model where we replace the uniformly distributed 
noise with Gaussian-distributed one 13511 



PUt)) 



1 



-«f(*)/V 



(6) 



where rj is the noise strength. In addition, in this version 
of the model particles with no neighbours are not sub- 
jected to angular noise. The repulsive interactions are 
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also noise-frec, so the term R{^i{t)) in Eq. ([T|) disap- 
pears. This kind of individual and collective behaviour, 
which implies some persistence in motion, is well suited 
for modelling biological systems and has been extensively 
used in literature jH, HJ, [36l - [39| . 

In both variations of the model the particle positions 
are updated by streaming along to new direction accord- 
ing to 



r,;(t + At) = r,{t) +V,{t + At) At 



(7) 



We use the standard Vicsek updating scheme (the for- 
ward updating rule, as defined in Refs. [H, HI])- 

In our calculations, we assume a constant propulsion 
speed vq = 1 and time step At = 1. We start the sim- 
ulation by randomly placing the particles inside the box 
and assigning random velocity directions. Total number 
of time steps is at least 2 x 10^ for all simulations. To 
set the required density we keep the box size constant 
{L = 500) and vary the number of particles in the inter- 
val from 250 to 150000. The statistics is collected in the 
steady state (after 50000 time steps) , and each character- 
istic of motion is calculated by averaging over 5 indepen- 
dent runs. As we mentioned above, the size of the zone 
of repulsion is set to unity in our simulations, Rr = 1. 
Moreover, we used the same radius of the alignment zone 
i?a = 5 in all runs. 



B. Motion statistics 

To analyse the collective behaviour of the model we 
perform a series of simulations changing the density, p, 
or 1], the noise strength. We characterise the collective 
motion with three different correlation/distribution func- 
tions. The velocity autocorrelation function is calculated 



V,(0)-V,(t) 



N \^ |V,(0)|-|V,(t)|/ 
The radial distribution function is defined by 

— \ 

U=l j^i 



N N 



(8) 



(9) 



where S is the Dirac delta function, (•) stands for the en- 
semble average. The two-point velocity correlation func- 
tion is calculated as 



TV N 



V^{t)■V,{t) 

N{N-1) \^^|V,(t)|.|V,(t)| 



v-1) 



(10) 



where i and j label particles separated by distance r = 
\rij\. With this definition, two particles with parallel (an- 
tiparallel) velocities give a correlation of -f-l (—1). The 
angular brackets denote the ensemble average. To char- 
acterise the swarming behaviour of the particles we also 



perform a cluster analysis. Cluster in our model is de- 
fined as a group of particles with a distance between 
neighbours smaller or equal to the radius of alignment 
zone Ra, i.e. particles interacting directly or via neigh- 
bouring agents are included into one cluster. This defini- 
tion is similar to one proposed in Ref . [l^ • We calculate 
the number of clusters and mean cluster size. 

We characterise the orientational ordering by the po- 
lar order parameter, which quantifies the alignment of 
the particle motion to the average instantaneous velocity 
vector 



ip{t) 



1 



Nvo 



N 



(11) 



This order parameter has been extensively used to de- 
scribe the orientational ordering in various systems of 
self-propelled particles [13, HH, [Ij, . It turns zero in 
the isotropic phase and assumes finite positive values in 
the ordered phase, which makes it easy to detect the tran- 
sition. However, at low densities it may be more diSicult 
to detect the transition due to relatively small number of 
particles constituting the system and large density fluc- 
tuations. To locate transition poin ts precisely we also 
calculated the Binder cumulant [41| 



Gr 



1 



(12) 



where (•)t stands for the time average and index L de- 
notes the value calculated in a system of size L. The 
most important property of the Binder cumulant is a 
very weak dependence on the system size so Gl takes a 
universal value at the critical point, which can be found 
as the intersection of all the curves Gl obtained at dif- 
ferent system sizes at fixed density. To detect the 
transition points in q ~ p plane precisely we plot three 
curves for different L at constant density and find the 
point where they cross each other. Then, we use those 
points to construct the phase diagram. 



III. RESULTS AND DISCUSSION 
A. Collective motion 

We performed simulations of a system of active par- 
ticles to fully map the region of the ordered collective 
behaviour on the p — ?/ plane. To illustrate the role of the 
both parameters we describe below the variation of the 
order along the iso-density and iso-?7 paths. 

First, we look at the behaviour of the system at con- 
stant density p ~ 0.04 while varying the intensity of 
angular noise 77. We expect the onset of ordered be- 
haviour at low noise levels and a disordered phase in the 
strong noise regime. Fig. [2] shows the alignment of parti- 
cles velocities at different noise strengths. At zero noise. 
Fig. ^a), interacting particles are aligned almost per- 
fectly. They can form massive clusters and orientation of 
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each individual completely coincides with the direction 
of movement of the group. With an increase of the noise 
amplitude the size of the correlated domain shrinks and 
at high levels of noise particles cannot form any stable 
clusters, Fig.[2l^b). We note that this kind of behaviour is 
an inalienable feature of models with metric interactions. 
Recent studies [12,111] on metric- free models showed that 
topological interactions cause significantly higher cohe- 
sion of the particles leading to formation of one or several 
large clusters with higher resistance to noise. 

Analysis of the velocity autocorrelation function in our 
model (Fig. [2Ka)) shows that at = motion of parti- 
cles is most persistent. The change of the direction of 
motion is realised only through collisions between differ- 
ent clusters. However, we see a fast decorrelation of the 
velocity in presence of strong angular noise. The qual- 
itative change from the slow decaying regime to a fast 
decay happens at 77 > 0.3 (Fig. [2Ia)). Plots of the radial 
distribution function. Fig. [31(b), and mean cluster size. 
Fig. m show that the particles tend to form clusters with 
the density of up to 3-4 times above the average. The 
cluster size can be as large as several hundred particle 
sizes (as set by Rr) and can include hundreds of individ- 
uals. The spatial range of the alignment extends as far as 
10 particle sizes at high noise {rj = 0.6, 0.8 in Fig. [2Ic)) 
but becomes much larger at the lower noise levels. On 
changing the noise, the decay also changes qualitatively: 
the spatial velocity correlation decays exponentially at 
high noise and with a power law at the lower noise levels. 

Figure U] presents the cluster statistics for the same 
density p ~ 0.04. The main plot shows the cluster size 
distribution on a log-log scale. We observe two quali- 
tatively different distributions: at low noise, rj = to 
0.3, the curves have a straight segment at large numbers, 
which indicates the power law decay of the distribution. 
The algebraic decay can be associated with a scale-free 
system where a formation of arbitrarily large clusters is 
possible. At 77 > 0.3, only exponential decay is observed 
and the formation of large clusters is suppressed. The 
plot in the inset shows the mean cluster size as a func- 
tion of the noise amplitude. Quite expectedly, the mean 
cluster size is decreasing with the noise amplification. 

We analyse now how the collective motion changes 
upon variation of the density at constant noise rj = 0.4, 
(Figs. El il). At low density (p = 0.004), the size of 
the correlated domain is minimal, particles form small 
clusters (Fig. [5]) that move in random directions. The 
clusters grow with the increase of the density and mo- 
tion becomes more persistent. Behaviour of the velocity 
autocorrelation function (Fig. [IKa)) at these conditions 
follows two distinct patterns: we see either fast decay 
within 10 to 30 time units (p = 0.004 or p = 0.5) or 
a very slow decay (all densities in between) such that 
there is no visible decay in the correlation function over 
hundred time steps. The former pattern is typical for 
stochastic disordered systems while the latter one corre- 
sponds to a very persistent and ordered motion. Plot of 




Figure 2: Simulation snapshots of a swarm in the Vicsek-type 
model (p = 0.04).|(a)|At zero noise. |(b)|At 77 = 0.6. 



the radial distribution function in Fig. |6l[b) shows that 
maximum local density of particles is observed within the 
zone of alignment but the increased density region spans 
to very large distances in the systems showing the per- 
sistent motion. We see a similar qualitative difference 
in the spatial velocity correlations for these systems, as 
shown in Fig. El^c): Two of them, p = 0.004 and p = 0.5, 
decay exponentially, while the other four show an alge- 
braic decay with C|| (r) oc r~'^ spanning over interval of 
more than hundred interparticle distances. 

In Fig. [71 we present the clustering statistics at vari- 
able density. Similar to the noise level series, we ob- 
serve two types of cluster size distributions: At low den- 
sities, p = 0.004 and 0.016, we see exponential distribu- 
tions, while the four other curves at p > 0.016 contain 
a straight segment, which corresponds to a power law 
decay. Moreover, at the three highest densities we also 
observe a peak at large particle numbers at m ~ 10000, 
which is close to the total number of particles in the sys- 
tem. We can conclude that particles travel in a single 
flock most of the time. The inset plot shows the mean 
cluster size as a function of the density. We note that at 
p < 0.024 no clusters are formed, while the cluster size 
starts growing rapidly at p > 0.1, i.e. there is a sharp 
transition into the "condensed" phase on increasing den- 
sity. We note that the cluster size distribution does not 
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Figure 3: Statistical properties of the Vicsek-type model at 
various noise magnitude (A'^ = 10000, p = 0.04). [(a)] Semi- 
log plot of velocity autocorrelation function C{t) over time 
t. (b) Radial distribution function g{r). 
correlation function C\\{r). 
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(c) Spatial velocity 



Figure 4: Cluster statistics for the Vicsek-type model at con- 
stant density (p = 0.04). The exponent p{m) oc m~'' for the 
straight segment: 77 = 0.1 - C « 0.77, = 0.2 - C f» 1.0, 
77 = 0.4 - C ~ 0.65. Inset: Average cluster size vs the noise 
strength. 



where authors found exponents ranging from C = 0.5 to 
C = 1.5 depending on the level of noise and the density. 
Although our simulation settings, the particle speed, and 
the size of the alignment zone differ from those in papers 
by Huepe et al. our values of C, lie in the same range 
[l8| . The differences in cluster size distributions between 
the Vicsek-type models and systems of rod-like particles 
[Slj might be related either to isotropic character of in- 
teraction or more long-range aligning interactions, thus 
leading to stronger cohesion in the former systems. 

From the results presented in this subsection we can 
conclude that motion in our model becomes persistent in 
a limited range of noise level and densities. Outside this 
region, the motion is predominantly random with short- 
range temporal and spatial correlations. The two types of 
collective dynamics are characterised by qualitatively dif- 
ferent statistical properties: exponentially decaying cor- 
relations in the disordered state and the algebraic decay 
and long-range correlations and scale-free behaviour in 
the ordered state. It is clearly seen that the persistence 
is a collective effect, which macroscopically corresponds 
to big groups of particles moving in the same direction 
and stabilising each other's direction of motion. We also 
see that in addition to orientational ordering a density 
separation similar to a vapour condensation into a liquid 
phase is taking place. The aligning interaction between 
particles acts in this case as an effective attraction, which 
suppresses the relative motion of the particles, thus lead- 
ing to particle aggregation. 



show the same exponent p{rn) oc to~^ as myxobactcrial 
systems or hard rods [sH ]. where an exponent Q of about 
0.9 has been found. Although the bulk of our measure- 
ments shows an exponent C = 0.9 1.1, there are some 
deviations at low densities, where the exponent becomes 
smaller (C « 0.66) and at high noise strength (C > 1.5). 
This can be also compared with the cluster size distri- 
butions obtained for the original Vicsek model [l^ 



B. Orientational ordering 

As one can see from the simulation snapshots, the ve- 
locity alignment in our model leads to global symme- 
try breaking, when most of the particles end up moving 
steadily in the same direction (Figs. HJa), [Sfb)). In this 
subsection we will analyse the transition in some more 
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Figure 5: Typical snapshots for the Vicsek-type model at 
various densities {r\ = 0.4). |(a)| p = 0.004. |(b)| p = 0.1. [(cj 



0.55. 



detail. 

The evolution of the orientational order parameter is 
illustrated by the time series in Fig. [5] for various levels 
of noise. At weak noise, 77 = — 0.2, the model demon- 
strates highly ordered motion with the order parameter 
fluctuating about the finite mean value. At the higher 
level of noise, 77 = 0.3, 0.4, we see a borderline behav- 
ior with large fluctuations, which are characteristic for 
the critical region, where the system fluctuates between 
the ordered and disordered states. Finally, at 77 = 0.6, 
the order parameter is zero within the statistical uncer- 



Figure 6: Statistical properties of the Vicsek-t ype model in 
presence of constant angular noise (77 



0.4). 



(a)| Semi -log 
(b) 



(c) Spatial velocity corre- 



plot of velocity autocorrelation fun ctio n G(£) over time t 
Radial distribution function g{r) 
lation function C\\(r). 



tainty. Fig. [9] shows the behaviour of the order param- 
eter for various values of the swarm density. Here we 
see a wildly fluctuating order parameter at p = 0.004 
and p = 0.04, then the finite values with some fluc- 
tuations at p = 0.1,0.2,0.3 and the absence of order at 
p = 0.55. In contrast with the dependence on the noise 
strength, where the ordering was gradually decreasing 
upon increase of the noise, the concentration behaviour 
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Figure 7: Cluster statistics for the Vicsek-type model in pres- 
ence of constant angular noise (77 = 0.4). The exponent 
p(m) oc m~'' for the straight segment: p = 0.04 - C ~ 0.66, 
p = 0.1 - C ~ 0.99, p = 0.15 - C ~ 1-1, P = 0.2 and p = 0.25 - 
( ~ 1.2. Inset: Average cluster size vs the density p. 
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Figure 8: Evolution the orientational order parameter in the 
Vicsek-type model at constant density (p = 0.04). The initial 
configuration at f = is disordered. 



is non-monotonous. The ordering is weak at small and 
high densities and high in between. The highest values 
of ip are reached at p = 0.1. This finding is in agreement 
with the observation made earlier (see Fig. 



5(c)): The 



spatial velocity correlation spans over smaller distances 
on increasing the particle number density. We can relate 
this to a larger concentration of orientational fluctuations 
at the higher densities as compared to the small system. 
The orientational fluctuations in the motion of the indi- 
viduals lead to particle overlaps and the temporary loss of 
alignment due to the repulsions, which leads to globally 
less ordered motion. 

It is clearly seen from the Figs. [8] and |9] that symmetry 
breaking into an orientationally ordered state is possible 
only after reaching certain levels of the noise strength 
and particle number density. We combined a large set 
of the order parameter data for various values of p and 
77 in Fig. [TUT c). The plot confirms the existence of the 
optimum intervals for both density and the magnitude 
of noise within which the global ordering (non-zero av- 



Figure 9: Evolution of the orientational order parameter in 
the Vicsek-type model at constant noise (77 = 0.3). The initial 
configuration at t = is disordered. 



erage i^) is maximal. The minimum density required for 
the ordering is quite small. In dense systems, the range 
of the ordered behaviour is widest at 77 = and is get- 
ting narrower upon increasing rj. We also note that the 
density, at which the maximum ordering is observed, is 
increasing from p w 0.1 at small 77 = to p « 0.18 at 
77 = 0.46. The stationary value of the order parameter in 
each configuration is determined by the competition be- 
tween the disordering noise and the aligning interactions. 
Therefore, the range of densities at which our model ex- 
hibits a transition to the ordered state depends on the 
level of noise. 

In Figs. [TUfa'l.fb). (p is plotted for various values of 
p and 77. At zero noise and low densities the order pa- 
rameter is close to 1, which corresponds to a completely 
aligned system. As the density and the noise increase the 
order parameter eventually decays to zero. We find that 
If varies continuously across the transition on varying the 
particle number density at the fixed noise level (see Fig. 
[S]and the map on Fig. Illf a)). Similarly, the transforma- 
tion that happens on increasing the noise level, when the 
density is fixed, is continuous with the order parameter 
decaying to zero (see Fig. [3] and the map on Fig. [TTTa)). 
in agreement with previous observation for the model 
with the same alignment and noise type [l^ [2^ . This 
conclusion is also confirmed by the behaviour of Binder 
cumulants across the transition. 

We should note that observed behaviour of the or- 
der parameter and correlation functions closely resembles 
those for Heisenbergferromagnet, as has been discussed 
previously [l^ [ij, [2^ [2^. An extensive discussion of 
the finite size scaling and critical ex pon ents in the Vic- 
sek model can be found in Rcfs. [23l. |29| . The change in 
the behaviour of C||(r), as seen in Figs. [3{c) andlHJc), 
can be associated with a symmetry breaking in the sys- 
tem and a transition into the aligned state. In the two- 
dimensional XY spin model, the decay of the two-point 
correlation function is exponential in the isotropic phase, 
while at the critic al p oint the decay becomes algebraic, 
Cii (r) oc [U l45l |. Note that generally the correla- 
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Figure 10: Orientational order parameter for the Vicsek-type 
model, (a) The iso-77 curves, (b) The iso-density curves, (c) 
The surface plot. 



tion function of the order parameter scales at the critical 
point as G(r) cx in d-dimensional systems j46j . 

In two dimensions, the simple spin- wave theory leads to a 
temperature dependent exponent r| = fcsr/(47rJ), where 
ks and T are the Boltzmann constant and absolute tem- 
perature and J is the spin coupling constant, while an 
inclusion of interaction between spin vortices yields an 
exponent t) = 1/4 [45l.[47|. In two-dimensional nematics, 
an algebraic decay of order parameter correlations have 
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Figure 11: Phase diagram for the Vicsek-type model, (a) 
Complete diagram in r; — p plane. The red and blue dots 
show the settings corresponding to series shown in Figs. [S] 
m and [51 [3 respectively, (b) Log-log plot showing transition 
region at low densities. 



been previously found across the whole nematic region, 
although with a temperature-dependent exponent [48| . 

In our system, the correlation function indeed decays 
exponentially in the symmetric (disordered) phase. In 
the non-symmetric (ordered) phase, it varies according to 



a power law C|| (r) oc 



where the exponent takes 



values from r| = 0.5 to t) « 0.97. At the transition point, 
however, the exponent is expected to satisfy the Fisher's 



scaling law: j/iy = 2 — t) 



46| . where 7 and v are the 



critical exponents for isothermal succeptibility and the 
fluctuation correlation radius. Indeed, in the limit of low 
concentrations, where the repulsions are not important, 
we have r| = 0.5, thus 2 — r| = 1.5, which is in agreement 
with the result for ^jv = 1.47 obtained previously for 
the standard 2D Vicsek model [2^. The other values of 
r| in our system are obtained inside the ordered region, 
far from the critical line, so the scaling relations are not 
expected to hold. 

In Fig. [12] we plotted the values of the exponent r| for 
two versions of the model (a) with uniformly distributed 
and (b) Gaussian noise (we consider it necessary to re- 
mind a reader that in this variation of the model particles 
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Figure 12: Behaviour of the scaling exponent of the velocity 
correlation function, C||(r) oc r~'^, for the Vicsek-type model 
with (a) uniformly distributed and (b) Gaussian noise for var- 
ious noise amplitudes 77 within the ordered phase. 
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Figure 13: Orientational order parameter for the Vicsek 
model with different radii of the repulsion zone Rr at vari- 
ous levels of the noise at p = 0.04 (a) and at variable density 
at 77 = 0.4 (b). 



with no neighbours as well as non-interacting agents are 
not subjected to the action of noise) [IHl- In contrast 
to the XY model; in the Vicsek-type system the value 
of T] seems to be insensitive to the "temperature" (the 
noise strength) and depends only on the density. We 
note that the lower values of the exponent are observed 
only at low densities, where the orientational phase tran- 
sition is accompanied by the gas-liquid-type transition. 
We can speculate that the crossover is reflecting the com- 
petition between the cluster size and the range of align- 
ment. The correlations are fully developed at densities 
p > 0.2, so the exponent remains constant at the higher 
densities. The higher value of the exponent in the Vicsek- 
type model as compared to the 2D XY is a consequence 
of the isotropy of the interaction, which represents the lo- 
cal mean field and is not bound to the lattice directions. 
Another important difference between the two models is 
that the active swarm with aligning interactions must 
be able do develop a true long-range order, i.e. the or- 
dered cluster can grow indefinitely in a large system as 
it follows from the nonlinear continuum theory 
The transition lines in the rj ~ p plane are best seen in 
Fig- fllTb). In the low-density region we observe the 
same power law 77 oc p'^ with k — 0.45 for both types 



of noise. This finding is in agreement with recent theo- 
retical and computational results for the original Vicsek 
model [21I, [22I, [2^. The phase diagram we observe in 
Figs. [111(a), (b) possesses one new important feature as 
compared to previously reported ones. We see the break- 
down of the global order at high densities. In our system, 
this breakdown is related to the misalignment caused by 
more frequent particle collisions as the system becomes 
more crowded. To confirm this hypothesis we have calcu- 
lated the collision frequency as a function of the density 
at various noise levels. Indeed, we found that the relative 
frequency of collisions continuously grows with the den- 
sity and reaches 0.8 at p = 0.6 (i.e. 80% of time the par- 
ticles are just bouncing back from each other). We also 
found that a removal of this mechanism or replacement 
of it by a more traditional excluded volume interaction 
removes the reentrant transition at the higher densities. 
In Fig. [13] we present results of simulations performed 
for various radii of the repulsion zone for the variation of 
our model with uniformly distributed noise. It is clearly 
seen that upon increase of Rr the transition point moves 
to the lower values of the noise amplitude (Fig. IT^ a)). 
At the same time, in the density scan (Fig. ITSf b)) the re- 
entrant transition into disordered phase happens earlier 
for the bigger radii of the repulsion zone and disappears 
at i?r = 0. We should note that the range of the or- 
dered behaviour on the phase diagram depends also on 
the propulsion speed and the time step and decrease of 
either of these variables leads to a shift of the transi- 
tion points to higher densities. Therefore, this destruc- 
tion of order seem to be a specific feature of the collision 
avoidance rule used in the Couzin's version of the model 
[251 . [33, S^l . We can imagine that in real life systems this 
kind of behaviour would be characteristic for swarms of 
agents with step-like motion such as high-density human 
crowds. 



IV. CONCLUSIONS 

We have studied the statistical properties of a system 
of self-propelled particles with aligning interactions and 
collision avoidance, representing a variation of the Vicsek 
model, across the dynamical order-disorder transition. 
We observed that in addition to the orientational transi- 
tion on increasing the particle number density reported 
earlier there exists a re-entrant transition into the disor- 
dered phase at high concentrations, which is related to 
the action of short-range repulsions destroying the local 
order. The transition in our system is continuous along 
both density and noise strength axes in agreement with 
previous observations of models with the same alignment 
and updating rules. We have demonstrated that orienta- 
tionally ordered phase is characterised by algebraic be- 
haviour of the two-point velocity correlation function and 
cluster size distribution. The critical exponent for the de- 
cay of the two-point velocity correlation function exhibits 
a crossover from t) = 0.5 to a density- independent value 
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r| w 1 on increasing density. 
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